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Spectral split in prompt supernova neutrino burst: Analytic three-flavor treatment 
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The prompt Ue burst from a core-collapse supernova (SN) is subject to both matter-induced flavor 
conversions and strong neutrino-neutrino refractive effects. For the lowest-mass progenitors, leading 
to 0-Ne-Mg core SNe, the matter density profile can be so steep that the usual MSW matter effects 
occur within the dense-neutrino region close to the neutrino sphere. In this case a "split" occurs 
in the emerging spectrum, i.e., the Ve flavor survival probability shows a step-like feature. We 
explain this feature analytically as a "MSW prepared spectral split." In a three-flavor treatment, 
the step-like feature actually consists of two narrowly spaced splits. They are determined by two 
combinations of flavor-lepton numbers that are conserved under collective oscillations. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

The neutrino flux streaming olT a collapsed supernova 
(SN) core is an intriguing astrophysical case where fla- 
vor transformations depend sensitively on some of the 
unknown elements of the Icptonic mixing matrix [l|, Q ■ 
Since the conversion probabilities depend also on the 
time-dependent matter profile, a high-statistics SN neu- 
trino observation may also reveaL for example, signatures 
for shock- wave propagation [1, [|, i, i, HSU- While 
galactic SNe are rare, various ongoing and future exper- 
imental programmes depend on large detectors that, be- 
sides their main purpose, are also sensitive to SN neutri- 
nos [To| . Therefore, understanding the flavor evolution 
of a SN neutrino signal remains of topical interest. 

The flavor transformation probabilities not only de- 
pend on the matter background, but also on the neutrino 
fluxes themselves: neutrino-neutrino interactions provide 



a nonlinear term in the equations of motion [ill, 12 1 
that causes collective flavor transformations [l3|, [ij, 
m, M, [13, M, M, S mi, [13. Only recently has 
it been fully appreciated that in the SN context these 
collective effectsgive rise to qualitatively ne w p henom- 

ena [H, Hi, [H, M M, M Mm, M, MMMMM- 

One peculiar aspect of the expected SN neutrino fluxes 
is the hierarchy F,,^ > Fa^ > F^^ = F^^ = = i^p^ so 
that there is an excess flux of Ve^e pairs over those of 
the other flavors. The nonlinear terms cause a collective 
transformation Vei^e — > VxVx^ where is a specific lin- 
ear combination of v^i and v^- The detailed dynamics 
of this transition is complicated and several important 
aspects are only numerically observed, not analytically 
understood. Still, the most crucial point is that the pair 
transformation V(,i>e ^ VxVx proceeds collectively much 
faster than ordinary pair annihilation, so we have to con- 
tend with a "speed-up phenomenon" [2l|, [l^l . The pair 
process does not violate flavor-lepton number. Being an 
instability in flavor space, it proceeds efficiently even for 
a very small mixing angle. 

A very different situation prevails in the interior of a 
SN core where the v,, distribution is determined by a large 



chemical potential that enhances the z^e density and sup- 
presses the Vf^ density relative to that of v^, and Vx so that 
collective pair conversions are not possible. Neutrino- 
neutrino interactions are strong, but their only impact 
is to synchronize the flavor oscillations ( "self-maintained 
coherence"). Significant flavor transformation here re- 
quires a violation of flavor-lepton number. However, the 
high density of ordinary matter suppresses the effective 
mixing angles between z/g and the other flavors. There- 
fore, in the interior of a SN core the individual flavor- 
lepton numbers are almost perfectly conserved [37|. 

Immediately after collapse, the SN emits a prompt v^^ 
burst that arises from the de-leptonization (neutroniza- 
tion) of the outer layers of the collapsed core. Once more 
we have a strongly enhanced and a suppressed v^. flux 
relative to the other flavors [s^. Once more, collective 
pair transformations are not possible: efficient flavor con- 
version requires a large violation of flavor-lepton number 
and is not possible if the ordinary matter density is large. 
At some distance from the neutrino sphere, the z/g flux 
encounters the usual MSW level crossings driven by the 
atmospheric neutrino mass difference Am^^jj^ (H cross- 




FIG. 1: Profile of the matter potential A and the effective 
neutrino-neutrino interaction potential /i for an O-Ne-Mg core 
collapse SN [H, 113, 53, nil . 
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FIG. 2: Mass eigenstate fractions Pa as well as the Ue survival probabilities far away from the star, numerically computed 
using the SN model of Fig. [T] and an initial flux of pure i/e . 



ing) and by the solar mass difference AtUq (L crossing), 
leading to the usual resonant transformations P, . 

An interesting new case is motivated by the insight 
that SNe with the lowest progenitor masses of 8-10 Mq, 
encompassing perhaps 30% of all cases, collapse before 
forming an iron core, the class of 0-Ne-Mg core collapse 
SNe [H, [13, \^ ■ In state-of-the-art numerical simu- 
lations these SNc explode even in a spherically symmet- 
ric treatment (convection plays no role), largely because 
their envelope mass is very small. By the same token, the 
matter density profile above the core is very steep even 
at the time of core bounce. In this case the H and L level 
crossings occur very close to the neutrino sphere and may 
well lie deeply within the collective neutrino region. This 
is illustrated in Fig. [1] where we show A(r) = \/2GFne{r) 
of an 0-Ne-Mg core progenitor star [H, |4l[ . We also 
show cjh = and wl = (AmQ/2i?) as hor- 

izontal lines, where the average is over the Fermi-Dirac 
spectrum of neutrino energies described below. The in- 
tersection of A(r) with these lines indicates the locations 
of the H and L level crossings. 

In Fig. [T] we also show the effective neutrino-neutrino 
interaction potential ^ = v^Gf-FI/^ (1 — cos6')cff, where 
9 is the angle between different neutrino trajectories and 
(. . .)eff stands for a suitable average. At large distances, 
H scales approximately as r~^. Collective neutrino effects 
driven by Arn^^^ are important for ^(r) > cjh and driven 
by Attiq for fj,{r) > wl- 

Duan et al. [U have recently shown that in this 
case the interplay of ordinary MSW conversions with 
collective oscillations leads to interesting effects. We 
start with a pure flux with a Fermi-Dirac spectrum 
{{E^^) = 11 MeV, degeneracy parameter 77 = 3), and nu- 
merically calculate the mass eigenstate fractions Pa and 



the h'e survival probabilities Pee far away from the star, 
as shown in Fig. (5] These plots are in qualitative a gree - 
ment with the corresponding curves in Fig. 2 of Ref. [3J] . 
However, our Pee is constructed as an incoherent sum of 
the mass fractions, thus representing the physical situa- 
tion far away from the star, where the oscillatory features 
seen in Duan ct al.'s Pee have disappeared. 

In inverted hierarchy, one observes that the neutrinos 
emerging from the star are in the V2 state at low energies 
and in the i^i state at high energies, with the transition 
taking place around i? « 12 MeV. This results in a step 
function in energy for Pee. 

In the normal hierarchy, the neutrinos emerging from 
the star are in state for i? > 17 MeV, in 1/2 state 
for 15 MeV < £: < 17 MeV, and in the 1^3 state for 
E ^ 15 MeV. The bump seen around 5 MeV is due to an 
abrupt change in the matter density profile used for the 
computation (see [s^l for details), and we do not address 
it here. The transition at i? « 15 MeV is rather sharp, 
however the one at « 17 MeV is not as abrupt. This 
results in a two-step function for Pee, with the step at 
E 17 MeV somewhat smoothened out. 

Broadly, this is an example of a "MSW prepared spec- 
tral split." In a two-flavor language, it is explained as 
follows. The strong neutrino-neutrino interactions lead 
to a synchronization of the neutrino oscillations. The 
flavor polarization vector of the ensemble begins at high 
density essentially aligned with the ly^. direction in flavor 
space. After passing the MSW region, the polarization 
vector emerges with a significant transverse component 
relative to the mass direction because the MSW tran- 
sition is not fully adiabatic. Subsequently this MSW- 
prepared initial condition is subject to collective effects 
only. As the effective neutrino-neutrino interaction be- 
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comes weaker, the modes above a certain energy -BspUt 
orient themselves along the mass direction, those with 
smaller energies in the opposite direction. This is pre- 
cisely the "neutrino only" case studied in Refs. [H, [2§| 
as a generic case for a spectral split, a case that requires 
one to prepare the polarization vector with a large trans- 
verse component. (For neutrinos plus antineutrinos, the 
MSW preparation is not necessary because the collective 
pair transformations alone engineer a split.) 

Our main goal here is to use the picture of a MSW- 
prepared spectral split to derive analytically the main 
features seen in Duan et al.'s numerical study [s^ . viz. 
the existence and the positions of the spectral splits. 

The numerical model in Fig. [1] shows that the MSW 
resonances are within the collective neutrino region, but 
not deeply inside. Therefore it is not a priori obvious 
if the "synchronized MSW preparation" and the subse- 
quent "split" can be clearly separated. For our discussion 
we therefore adopt a more schematic model. We artifi- 
cially increase the neutrino-neutrino interaction strength 
(raise the ^-profile in Fig. [1]) such that the MSW region 
and the spectral split regions are clearly separate. In this 
framework we first calculate the MSW preparation ana- 
lytically and then study a three-flavor treatment of the 
spectral split, based on the machinery recently developed 
by two of us [1^ . We find that our analytic treatment 
reproduces the numerical results surprisingly well. It also 
reproduces all relevant features of the realistic case, i.e. 
with the /i-profile of Fig. [U justifying our simplifying as- 
sumptions and verifying our general interpretation. 

We first set up in Sec. [11] our schematic SN model 
that captures the features relevant for our treatment and 
continue in Sec. IIIII with the equations of motion. In 
Sec. IIVI we derive analytically the MSW-prepared three- 
flavor state that serves as input for our three-flavor spec- 
tral split study in Sec. El We conclude in Sec. |VT1 

II. SIMPLIFIED SUPERNOVA MODEL 

A. Realistic scenario 

We take the electron density proflle rip, obtained from 
numerical studies of 0-Ne-Mg core SNe j39l . [40l . [4l| , pro- 
viding A(r) as shown in our Fig. [TJ The neutrino lumi- 
nosity is assumed to be L^^ ~ 10^^ erg s"-'^ with a Fermi- 
Dirac spectrum with the average energy (E^^) = 11 MeV 
and a degeneracy parameter 77 = 3, implying a temper- 
ature T = 2.76 MeV. The neutrino sphere is taken at 
the radius i? = 60 km, implying an effective neutrino- 
neutrino interaction strength at large distances of 

Mr) = Mo(^)' (1) 

with fiQ = 8.6 X 10^^ kni^^. This is the profile shown in 
Fig. [TJ Based on this model we have solved the equations 
of motion numerically and found the electron neutrino 



survival probability shown in Fig. [2l in agreement with 
the results of Duan et al. [11]. 

In this realistic situation, the decrease of the effective 
neutrino-neutrino interaction with radius is such that the 
spectral split is essentially adiabatic. Quantitatively, the 
length scale £^ = \d\n fi{r)/dr\~^ is large enough to sat- 
isfy the adiabaticity condition (see Sec. V of [2^) during 
the spectral split. The very development of a step-like 
feature in P^e in Fig. [51 is an expression of the adiabatic- 
ity. The "sharpness" of the step in Pee is a measure of 
the degree of adiabaticity, an extremely slowly decreas- 
ing neutrino-neutrino interaction corresponds to a perfect 
step function in P^e 



B. Analytical treatment 

Our analytic treatment is based on a schematic repre- 
sentation of the essentials of this realistic case. We pic- 
ture that the synchronized MSW effect occurs first and 
factorizes from the collective oscillations that lead to the 
spectral split. To this end, we consider the limiting case 
where fiQ in eq. ([1]) is arbitrarily large. This makes the H 
and L level crossings more synchronized, and pushes the 
regions where the spectral splits occur, i.e. the regions 
where /i(r) becomes comparable to wh or wl, to much 
larger radii. This also renders the spectral splits even 
more adiabatic ^, making the steps in Pee sharper. 

It is also important that the MSW transition is not 
perfectly adiabatic. In our schematic model, we use the 
power-law profile for the matter potential, 

A(r) = Ao(^)" (2) 

with Ao = 10^ km~\ ro = 900 km and a 50. This ap- 
proximates reasonably the numerical A(r) profile of Fig.[T] 
in the neighborhood of the H and L crossings. This allows 
an analytic estimate of the level crossing probabilities at 
the two resonances. 

Thus, we obtain our analytic results in the limit of 
/io — *■ 00 and a power law profile for A(r), where (i) the 
MSW transitions are perfectly synchronized but semia- 
diabatic, (ii) the spectral split regions are well separated 
from the MSW region, and (iii) the spectral splits are per- 
fectly adiabatic, making the steps in P^e infinitely sharp. 
Apart from these minor changes, the physical pictures 
with the realistic profile and our schematic profile are 
identical. 



^ Since /i(r) is a power law, the length scale increases with 
increasing r. As a result, larger radius for spectral split implies 
larger and more adiabaticity 



4 



C. Numerical treatment 

We also study the spectral splits numerically, assuming 
the analytic MSW-prepared spectra as input. For such a 
numerical illustration we employ 



I^lq = 10^^ km' 



(3) 



much larger than the realistic value shown in Fig. [TJ 

The analytic results match the numerical ones in all 
details, and reproduce the main features of the realistic 
situation, with a minor difference that the spectral splits 
appear sharper. This is because the net effect of our ap- 
proximations is only to make the MSW resonances more 
synchronized, and the spectral splits more adiabatic. 



III. EQUATIONS OF MOTION 

A. Matrices of density 

Mixed neutrinos are described by matrices of den- 
sity Qp for each momentum. The diagonal entries are 
the usual occupation numbers whereas the off-diagonal 
terms encode phase information. The equations of mo- 
tion (EOMs) are 



idtQp = [Hp,£»p] , 
where the Hamiltonian is 11211 



(4) 



Hp = f7p+V+x/2GF J ^ {gci - Pq) (1-Vq-Vp) , (5) 

Vp being the velocity. The matrix of vacuum oscillation 
frequencies is fip = diag(m^, to|, m|)/2|p| in the mass 
basis. The matter effect is represented, in the weak in- 
teraction basis, by V = \/2 Gp^^e diag(l, 0, 0). While in 
general there is a second-order difference between the z^^ 
and Vr refractive index [43j that can be important for 
collective neutrino oscillations [s^, for the low matter 
densities relevant in our case this "mu-tau matter term" 
is irrelevant. 

The factor (1 — Vq • Vp) in Hp implies "multi-angle 
effects" for neutrinos moving on different trajectories 
[21I . [2^ . [23 |. However, for realistic SN conditions the 
modifications are small, allowing for a single-angle ap- 
proximation [2^[30|. In the strongly synchronized regime 
this is not surprising: the strong neutrino-neutrino inter- 
action causes self-maintained coherence not only between 
different energy modes, but also between different angu- 
lar modes. It has not been explained why the single-angle 
approximation remains good even when the neutrino- 
neutrino interaction becomes weak, although numerically 
this is observed to be the case [H, [M, ■ 

We are studying the spatial evolution of the neutrino 
fluxes in a quasi-stationary situation. Therefore, the ma- 
trices Qp do not depend on time explicitly so that the to- 
tal time derivative in the EOMs reduces to the Liouville 



term involving only spatial derivatives. (See Ref. [4J| for 
a recent comprehensive discussion of the role of the Li- 
ouville term in Boltzmann collision equations involving 
oscillating neutrinos.) Moreover, we consider a spheri- 
cally symmetric system so that the only spatial variable 
is the radial coordinate r. In the single-angle approxima- 
tion we finally need to study the simple EOMs, 

idrQoj = [H„, g^] , (6) 
where we now classify different modes by the variable 



Lu — . 

2E 

The single-mode Hamiltonians are 



H„ 



+ X{r)L + fi{r) g . 



(7) 



(8) 



Here we have introduced A(r) = ^/2Gpnf,{r) and L = 
diag(l,0, 0) in the weak interaction basis. The matrix 
of the total density g = J duj guj is normalized such that 
at the neutrino sphere it is diag(l,0,0) in the weak 
interaction basis. It is conserved except for oscillation 
effects, i.e., the physical neutrino density has been ab- 
sorbed in the coefficient that measures a suitable 
angular average of the neutrino-neutrino interaction en- 
ergy. The radial variation is /x(r) oc r^^ where a factor 
r^^ comes from the geometric flux dilution, another ap- 
proximate factor from the fact that neutrinos become 
more collinear with distance from the source. 



B. Mixing parameters 

Since Hi^ appears in a commutator we may arbitrar- 
ily add terms proportional to the 3x3 unit matrix. It 
will prove convenient to express the matrix of vacuum 
oscillation frequencies in the form 



il^ ^ LU diag (— ^a, 



ia,±l) 



where the mass hierarchy parameter is 

Ami, 1 
30 ■ 



(9) 



(10) 



A positive sign in the third component of signifies 
the normal mass hierarchy, a negative sign the inverted 
hierarchy. For the mass differences themselves we use [45j 



Ato^ = 2.4 X 10~3 eV^ . 



8 X lO^'^ eV^ 
For the mixing angles we use 

012 = 0.6, 

013 = 0.1. 



(11) 



(12) 



Of course, for ^13 only upper limits exist. In this con- 
text the CP phase Sqp can be ignored, since it does not 
influence the relevant probabilities for equal and i^^ 
fluxes HEl. 
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C. Bloch vectors 

111 the two-flavor context it is well known that the ma- 
trices of density can be expressed in terms of Bloch vec- 
tors, leading to EOMs that resemble the precession of a 
gyroscope around an external force field. This picture 
helps to recognize properties of the EOMs that are dif- 
ficult to fathom in the commutator form of the EOMs. 
Therefore, we follow a recent paper by two of us [s^ 
and note that every Hermitean 3x3 matrix X can be ex- 
pressed in the form 



X = ixo + ix.A, 



(13) 



where Xq — Tr(X), X is an eight-dimensional Bloch vec- 
tor, and A is a vector of the Gell-Mann matrices. Note 
that A3 = diag(l,-l,0) and Ag = diag(l, 1, -2)/a/3. 
The Bloch vector for the single-mode Hamiltonian is 
whereas the one for g^^ is the polarization vector P;^ . 
Then the single-mode EOMs are 



(14) 



where the cross product is understood in the SU(3) sense: 
(A X B),; = fijkAjBk where i, j,k^ 1, . . . , 8 and fijk arc 
the structure constants of the SU(3) Lie algebra. 

With the global polarization vector P = j duj P^j and 
ignoring the ordinary matter term, we write the single- 
mode Hamiltonian in the form 

= Lo {Bn + aBi,) + ^lP . (15) 

In the mass basis, the "magnetic field" components are 
2 

BH==-;^e8 and Bl = -63 , (16) 

representing the "atmospheric" and "solar mass direc- 
tions," respectively. 

In the absence of ordinary matter, the EOM for the 
global polarization vector is 



P = (Bh + aBL) X M, 



(17) 



where the "magnetic moment" of the system is M = 
J duj uj Pi^ . In the mass basis this is 



2 



a 63 I X M . 



(18) 



The vector on the r.h.s. is orthogonal to both 63 and eg. 
The reason is that the matrices A3 and Ag commute or, 
in other words, that fass = for a = 1, ... ,8 and the 
same for all permutations. As a consequence, the vector 
P has no 63 or eg component so that P3 = and Pg = 0. 
In a general basis this implies 



dt(P-BH)=0 and dt(P-BL) = 0. 



(19) 



This is the equivalent of "flavor-lcpton number conserva- 
tion" dt(P • B) = in the two-flavor context [H,!!!,!!^. 



In other words, in the three-flavor context we have two 
flavor-lepton numbers that are separately conserved. In 
the mass basis one concludes that 



^3 

Ps 



911 
911 



Q22 , 
- Q22 



2^3 



(20) 



are conserved. 



IV. SYNCHRONIZED MSW EFFECT 

As a first step in our analytic study we consider the 
matter-induced conversion of the initial fiux as it 
passes the H and L level crossings, assuming that both 
lie deeply in the region where the neutrino-neutrino in- 
teraction is strong. As a result, the flavor oscillations arc 
synchronized, meaning that all 9uj stay pinned to each 
other. In other words, it is enough to study the evolution 
of the matrix of the total density 9. It obeys the EOM 



idr9{r) = [n + \{r)L,g{r)] 



where in the mass basis 



i7 = (jjo diag (— ^a, +\a, ±l) 



(21) 



(22) 



&nduJo = {^ml,j2E). 

For the Fcrmi-Dirac spectrum described in Sec. [TT] we 
find ujQ = 0.710 km^^. We will also use the notation 



2E 

,2 



a>o = 0.710 km"^ , 



=ac^o- 0.024 km- ^ (23) 



The matrix of vacuum oscillation frequencies thus can 
also be written as = diag (— ^wl, i^^^n) • 

The assumed perfect synchronization implies that we 
can treat this system as an equivalent system with a sin- 
gle energy^ or rather with two fixed vacuum frequencies 
tJH and determining Jl. It is most useful to study its 
evolution in the basis of instantaneous propagation eigen- 
states where we denote the total matrix of density as 9. 
Since the only effect of the neutrino- neutrino interactions 
is to synchronize the oscillations and to reduce the sys- 
tem to an equivalent single-energy case, the propagation 
eigenstates are defined by the ordinary matter term for 
a monochromatic neutrino beam, wheras the neutrino- 
neutrino interaction plays no further role. The propaga- 
tion basis coincides with the weak interaction basis when 



^ An MSW transition in the presence of a dense neutrino gas was 
first treated in this way in the context of early-universe neu- 
trino oscillations by Wong [Tsll and by Abazajian, Beacom and 
Bell [l9l |. They used the terminology "collective MSW-like trans- 
formation" and "synchronized MSW effect," respectively. 
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Normal hierarchy 



Inverted hierarchy 




10 10' 



10 10^ 10"^ lo"^ 



Wo/A(r) 



FIG. 3: Evolution of the diagonal elements of the matrix of densities q^i for a mode with lu — 0.42 km ^, represented in the 
basis of instantaneous propagation eigenstates. Due to the synchronization, all modes behave similarly. 



TABLE I: Final occupations of the mass eigenstates after 
passing both MSW level crossings as described in the text. 





Normal hierarchy 


Inverted hierarchy 


gii 


Ph Pl 0.12 


Pl 0.31 


922 


Ph (1 - Pl) 0.26 


(1 - Pl) 0.69 


Q33 


(1 - Ph) 0.62 






the matter density is large so that our initial state is 
g ~ diag(0, 0, 1) in normal hierarchy, and diag(0, 1, 0) in 
inverted hierarchy. If the subsequent evolution were per- 
fectly adiabatic, the system would remain in this state so 
that the j/g survival probability after the H and L cross- 
ings would be given by well known results 

We assume indeed that the evolution is adiabatic, ex- 
cept near the H and L level crossings where in our sys- 
tem the jumping probabilities Ph and Pl need not be 
small. The neutrino mass-gap hierarchy ensures that the 
two crossings factorize with good approximation. In the 
normal hierarchy, the system encounters both crossings 
and the final occupation of the propagation eigenstates 
is given by the products of probabilities shown in Ta- 
ble IT] In the inverted hierarchy, only the L crossing is 
encountered because the H level crossing is now in the 
antineutrino sector that is irrelevant in our case. Again, 
the final mass-state occupations are shown in Table [D 

The jumping probability for incomplete adiabaticity is 
given to a good approximation by the so-called double- 
exponential formula [i^, H^, [S^l 



Ph 



cxp(27rPHWH cos^ Ois) — 1 



where the scale height is 

d\nX{r) 



dr 



(24) 



(25) 



It has to be evaluated at the point of maximum violation 
of adiabaticity [5l|, [H, [11] , given by = A(rH) The 



assumed power-law profile of Eq. Q and our choices for 
wh and imply rn = 1089 km, Ph = 21.8 km, and 
Ph = 0.38. Analogous results pertain to Pl with H ^ L 
everywhere and the substitution ^13 612 . We find tl = 
1166 km, Pl = 23.3 km, and Pl = 0.31. These numerical 
results for Ph and Pl imply the numerical results for the 
final occupations shown in Table |T1 

In the steep density profile used here, the H and L 
resonances look spatially very close (Fig. [T]) so that one 
may worry if they indeed factorize in the usual way. We 
stress that the two resonances do not overlap, but one 
may still worry about possible interference effects. How- 
ever, we can compare the analytic results with a numeri- 
cal three-flavor evolution, using the same power-law A(r) 
profile. In Fig. [3] we show the evolution of the diagonal 
elements of g^j, for a mode with uj = 0.42 km"-'^, as a 
function of Uo/Xir). Due to the synchronization condi- 
tion, all frequency modes behave the same. On this scale, 
the H crossing is at ujo/X{r) = 1 and the L crossing at 
cjo/A(r) = a^^ = 30. The agreement between the nu- 
merical end states and the analytically predicted ones is 
striking. 



V. SPECTRAL SPLIT 

After the system has passed the two MSW level cross- 
ings, the ordinary matter density quickly becomes neg- 
ligible, whereas by assumption the neutrino-neutrino ef- 
fects arc still strong. The subsequent evolution to the 
point where the neutrino-neutrino interaction becomes 
negligible will then produce spectral splits in the same 
way as described in Refs. [2l,|2^. We can follow the pre- 
vious two-flavor treatment almost step by step because 
the present three-flavor system is simplified by the mass- 
gap hierarchy a w 1/30 ^ 1. While the two conserved 
flavor-lepton numbers present in the three flavor case lead 
to two spectral splits, these will occur in sequence and 
their dynamics factorizes in practice. 

The first split to develop is driven by the atmospheric 




FIG. 5: Evolution of the 3 and 8 components of Pi^ for our box-spectrum example. 
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mass difference and thus can be called the H split. As in 
Refs. m, [2§| we go to a rotating frame, at first rotating 
around the Bh direction. The single-mode Hamiltonians 
in this co-rotating frame are 



{uj -UjIi)Bu+ ^iP : 



(26) 



neglecting for now the much smaller term owBl- This 
is justified because, when fi ^ cu (and thus /i 3> aw), 
the ensemble of neutrinos is in a regime where we expect 
spectral splitting along eg and synchronized oscillations 
along 63. This factorization has been explicitly shown in 
[3^. Flavor conversion is thus driven primarily by Bh, 
while Bl gives sub-leading corrections due to the syn- 
chronized oscillations. Similarly, when /z ~ auj, flavor 
conversion proceeds efficiently via a spectral split along 
63 and is driven by Bl, while Bh drives vacuum oscilla- 
tions along eg. 

Now, as ^ adiabatically goes to zero, the co-rotation 
frequency c^h approaches the final split frequency cjfj ^^'^ 
the modes with uj > cu^ will orient themselves along Bh, 
those with lu < ujf^ in the — Bh direction. The value of 
wfj is fixed by the conservation of Pg. Since the evolution 
associated with Bh has saturated, we can next go into a 
frame rotating around Bl where 



(cj-cj£)BL + /iP 



(27) 



and repeat the analogous argument. 

To illustrate the dynamics of the split in a form similar 
to Refs. ,28i . .29i | wc consider an explicit example with an 
initial "box spectrum" at high density of the form 



e(c.) = {(2-o; 



-1 



for < a; < 2^0 , 
otherwise . 



At high densities Pee{^) coincides with p33(w) in normal 
hierarchy and with p22{^) in inverted hierarchy. After 
the MSW crossings the spectrum is still of box shape 
because of the assumed strong neutrino-neutrino inter- 
action, but now has the qu, Q22, and ^33 components 
shown in Table [J Note that after the MSW transitions 
we neglect ordinary matter so that the propagation eigen- 
states are identical with the mass eigenstates and g = g. 



A. Normal hierarchy 

At first we study the evolution caused by the neutrino- 
neutrino interactions numerically. To this end we use the 
usual /i(r) oc profile of Eq. ^ with the coefficient /ig 
of Eq. ([3]) . In this way the evolution is strongly but not 
perfectly adiabatic. Of course, the analytic results, based 
on the conservation of flavor-lepton number, apply in the 
perfectly adiabatic limit which here requires /io — *■ 00. 

In Fig. d] (left column) we show the evolution in the 
mass basis of the diagonal elements of p^j for 50 modes 
as a function of loq/ p.- All modes start with the same ini- 
tial condition prepared by the MSW transitions. Around 
ujq/p = 1 one recognizes a first H split that affects all 



components. Around ujo/p = one recognizes a sec- 
ond split, the L split, that affects only pu and p22- Dur- 
ing the H split, the modes with w < wfj tend to P33 0, 
while those with lu > approach P33 — > 1. Then pn 
and P22 modes with uj > go to 0, as implied by the 
conservation of the trace of p. While the modes with 
LU < rise towards higher values of pu and P22, they 
encounter the L split at a frequency ujf^ < w^j. At this 
split, for UJ < tof^, the pn approach 1 and the P22 ap- 
proach 0, and vice versa for < w < wfj. As a result 
of imperfect adiabaticity some modes do not reach these 
extreme values, but get frozen earlier. 

In Fig. [5] (left column) we show the same case in terms 
of the mass-basis P^i^z and P^^g components. We observe 
that in P^^.g only the H split operates whereas in P^,3 the 
H and L splits operate in sequence. 

The situation can be visualized in terms of the e3-eg 
triangle diagram [3^ shown in Fig. [S] Each point in the 
interior and on the boundary of the triangle represents 
the projection of the polarization vector in the e3-eg 
plane. Neutrinos from the Ve burst start in the state ~ 
, where by "tilde" , we represent the instantaneous mass 
eigenstates. The H crossing shifts the neutrino state from 
the vertex towards the i>2 state, but only partially, due 
to the semiadiabatic nature of the transition. After that 
crossing, all neutrinos find themselves at the point A' 
inside the triangle. The L crossing further transports the 
state along a line parallel to the i'2~i'i edge towards ui^ 
again only partly due to the semiadiabaticity. Before the 
split, all the neutrinos arc thus at a point A in the interior 
of the triangle. 

The H split takes the oj > ojf^ modes towards the j/3 



(28) state (Paj,3 = 0, Pi^^g = — 2/\/3) and the modes u < 



towards some combination of I'l and i/2, while conserving 
the total P3 and Pg. Since a ^ 1, the H and L splits are 
well separated and the high-w modes reach the i^s vertex, 
i.e. the H split saturates, before the L split begins. The 
\ow-uj modes propagating towards the P^^.g = 1/V3 line 
encounter the L split that tends to take the lu > ojf^ 
modes towards D2 (Pw,3 = — l,Pw,8 = 1/V3) and the 
LU < Luf^ modes towards Di {Pi^,3 ~ l,Pw,8 = I/a/S)- In 
the adiabatic limit, given sufficient time to propagate, 
the H and L splits result in all neutrinos reaching one of 
the three vertices of the e3-eg triangle. 

Using the conservation law for P3 and Pg of Eq. (PO]) 
one can evaluate the split frequencies wfj and a;£. For 
LU < wfj we have P^^^^g l/\/3, while for cu > wfj they 
reach —2/ a/3. In the limit of perfect adiabaticity, the 
conservation of Pg implies 



71^ "7! 



(20-0 



(29) 



where P° g is the common value of P^j^s before the split 
be gins. In our example, g = —0.50, leading to lu^ = 
0.76 ujo- 

When the H split saturates, all modes with lu > cj^j 
have P^j.g ~ —2/^/3, and hence P(^^3 = due to the 



9 



Normal hierarchy 

ec 



Inverted hierarchy 



V2 (-1,1 /V3) 




Yi (1,1/V3) 



V3 (0,2/V3) 



V2 (-1,1/V3)_ 


Vi (1,1/V3) 








T 63 
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FIG. 6: Projection of the polarization vectors P^i on the es-eg plane for our box-example. The vertices of the triangle represent 
pure (instantaneous) mass eigenstates. After both MSW transitions, the system is at the point A in the interior of the triangle. 
(See the text for details.) 
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FIG. 7: Diagonal elements of Quj. Thin line: initial box spectrum. Thick line: final numerical spectrum. Dotted line: final 
spectrum in the adiabatic limit. 
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FIG. 9: The Ve component at a large distance so that the different u modes have kinematically decohered. Convention for 
lines is same as in Fig. [T] 
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FIG. 10: Initial (thin) and final (thick) spectrum for a Fermi-Dirac distribution with the parameters described in Sec. Ull The 
numerical final spectrum is for our toy-model SN where the MSW crossings and spectral-split region are far separated. Dotted 
curves represent the survival probability Pee for electron neutrinos. 
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conservation of the norm of P. These modes have reached 
the bottom vertex of the e3-eg triangle and hence cannot 
spht further due to the L spht. On the other hand, for 
modes with lo < wfj a second split in P^^^s happens. These 
modes approach P^^ ^ = +1 for oj < ujf^ and P^^^ = —1 for 
Lu > Lof^. Applying the conservation law for P3 gives us 



(30) 



In our example P° 3 — —0.14 so that ujf^ = 0.24 wq- 

We show in Figs. [7] and [8] the mass-basis spectra of the 
diagonal elements of g and of P3 and Pg, . Thin lines are 
the MSW-prepared initial spectra. Thick lines show the 
numerical end states, corresponding to the split diagrams 
of Figs. [4] and m Dotted lines show the adiabatic limit- 
ing behavior based on the lepton-number conservation 
laws. Once more the agreement is striking. Imperfect 
adiabaticity leads to a smoothening of the splits which 
otherwise are sharp spectral steps. 

Finally, in Fig. [5] we show the spectrum after the 
two splits. The solid curve is the numerical result, where 
the vacuum oscillations between the SN and the observer 
have been averaged (kinematical decoherence between 
different uj modes), i.e.. 



of the box and thus is not visible. The conservation law 
for P3 and using in our case 3 = —0.38, one obtains 
t^L = 0.62^0. For the electron flavor we predict for the 
final spectrum 



Pe 



cos^ 9i2 for UJ < ujf^, 
sin^ 6*12 for oj? < iv . 



(33) 



9e 



UliQll 



UI2Q22 + C/e3P33 



(31) 



in agreement with the numerical result. 

Note that in the inverted hierarchy, there is only one 
split. This is because P^^^ is already at an extreme value 
before the split can begin. In general, there are two splits 
if the neutrino state before the split is in the interior of 
the 63-88 triangle, one split if it is along one of the edges 
of the triangle (as in this case), and no split occurs if the 
neutrino state is at any of the three vertices (See Fig. [6]). 



C. Fermi-Dirac spectrum 

It is straightforward to extend these arguments to a 
general spectrum, e.g. a Fermi-Dirac spectrum /{lo). The 
conservation laws imply for the normal hierarchy 

V3P" - / diuficu) -2 / dcuf{Lu), 



The dotted curve is our analytic result in the adiabatic 
limit. This result is easily explained if we observe that 
in Eq. (|3ip in the case of maximal 23 mixing one has 



^el 

sin^ 



= cos fi3 cos" (^12, 
hs- Therefore, 



U^2 = cos^ sin^ 



712 , and = 



for UJ < LJ 
for 

for UJ? 



L ' 



(32) 



< UJ . 



Again the agreement between the analytic and numerical 
results is very good. 



B. Inverted hierarchy 

For the inverted hierarchy we show the analogous in- 
formation in the right-handed columns of Figs. ISHHl The 
initial state here is !>2- The nonadiabatic L crossing takes 
the neutrino states partly towards i^i . After the L cross- 
ing and before the split, the neutrino state for all modes 
is along the i'i~i'2 edge, at A as shown in Fig. [S] (right 
column), where P^^s = I/a/S- Since all neutrinos al- 
ready are in one of the extreme values of Pi^.s-, the H 
split is inoperational. This corresponds to P33 remain- 
ing in its MSW-prepared initial value of 0. The L split 



-1 for w < Wl and pn for uj > uj- 



takes pii - 

and vice versa for p22- In the inverted hierarchy we have 
an effective two-flavor case in the z/1-1/2 subsector. This 
is a consequence of the MSW-prepared initial condition. 
Initially P^ g, ~ ^/ V^- Applying now the conservation of 
Pg we obtain wfj — 2ujo, i.e., the split occurs at the edge 



P^. 



duj f{uj) 



dujfiuj). (34) 



On the other hand, for the inverted hierarchy we find 



pO _ 



dujf{ui) — / du! f{uj) 



and ujf^. 



(35) 



Note 



These relations allow us to calculate cj^j 
that these results are exact only in the limit of an infinite 
mass-gap hierarchy, i.e., for a 0, where the H and 
L splits perfectly factorize. Once the split frequencies 
have been found, the corresponding energies are Ef^ = 
AmlJ^u;!, and El = Amy2a;£. 

For our schematic SN model where the MSW level 
crossings and the spectral split region are widely sepa- 
rated, we find El^ = 11.9 MeV and El = 16.9 MeV in 
the normal hierarchy, and El = 12.7 MeV in the inverted 
hierarchy. We show the initial and final spectra for the 
Fermi-Dirac case in Fig. 1101 Once more, we have coarse- 
grained over neighboring modes, representing the effect 
of kinematical decoherence. 



VI. CONCLUSIONS 

We have studied the three-flavor evolution of a i^e burst 
that first undergoes two MSW level crossings driven by 

An ~ 



J and Attiq, respectively, and then undergoes spec- 
tral splits by the adiabatically decreasing strength of 
the neutrino-neutrino interaction. This case study of an 
MSW-prepared spectral split serves as a proxy for the 
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recent numerical study of the prompt v^^ burst in an O- 
Ne-Mg core collapse SN. Here, the matter density profile 
is so steep that the sequence between MSW crossings 
and collective neutrino oscillations is reversed from what 
would be expected in a traditional iron-core SN. 

First we have analytically estimated the population of 
the propagation eigenstates after the MSW transforma- 
tions. Due to the sharply falling matter density at the 
edge of the core of the star, the MSW transitions arc not 
completely adiabatic, which helps in satisfying a precon- 
dition for the spectral splits to take place. We have used 
the well-known analytic double-exponential formula for 
calculating the jump probabilities. 

We have studied numerically the subsequent spectral 
splits. We have analytically determined the split frequen- 
cies on the basis of two conserved flavor-lepton number 
combinations that supersede the single conservation law 
encountered in a two-flavor situation. The neutrino mass- 
gap hierarchy allows for a factorization of an H split and 
an L split, similar to the factorization of the MSW effect 
into an H and an L crossing. The dynamics of the split 
evolution is clearly seen to be a two-step process. 

The dynamics of the two spectral splits can be under- 
stood in terms of the motion of the neutrino state in the 
63-63 triangle diagram, which can explain many of the 
features of neutrino evolution qualitatively. The number 
of possible splits can be deduced by the location of the 
neutrino state inside the triangle. We have also shown 
how the positions of the splits can be calculated accu- 
rately given the initial neutrino spectra, and calculated 
the i/e survival probability analytically, that matches the 
numerical computations. 

Our analytic treatment accounts very nicely for the 
numerical findings of Duan et al. [s^l • In their case the 
MSW conversion and the spectral splits are spatially very 
close so that it is not a priori obvious that our schematic 
model would be a good representation. In our treatment 
we have enforced a clear separation between the MSW re- 
gion and the spectral-split region by assuming the limit of 



large neutrino-neutrino interactions, a limit that also en- 
sures that the MSW transition is perfectly synchronized 
and that the spectral split is perfectly adiabatic. A pos- 
teriori, however, our interpretation as a MSW prepared 
spectral split appears nicely justified and quantitatively 
appropriate. Our treatment is also a useful application of 
the three-flavor oscillation machinery developed by two 
of us recently [sll . 

It appears that the impact of collective neutrino oscil- 
lations on the propagation of the prompt v^, burst is con- 
ceptually and quantitatively well under control. Char- 
acteristic signatures of these flavor transitions in large 
under grou nd detectors have also been recently investi- 
gated [54| . What remains is to observe these features in 
the neutrino signal of the next galactic SN. 

Note Added 

After our manuscript was completed, a paper by Duan, 
Fuller and Qian appeared that treats three-flavor split 
phenomena 
work. 
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